Non-abelian statistics from an abelian model 
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It is well known that the abelian Z2 anyonic model (toric code) can be realized on a highly 
entangled two-dimensional spin lattice, where the anyons are quasiparticles located at the endpoints 
of string-like concatenations of Pauli operators. Here we show that the same entangled states of the 
same lattice are capable of supporting the non-abelian Ising model, where the concatenated operators 
are elements of the Clifford group. The Ising anyons are shown to be essentially superpositions of the 
abelian toric code anyons, reproducing the required fusion, braiding and statistical properties. We 
propose a string framing and ancillary qubits to implement the non-trivial chirality of this model. 

PACS numbers: 05.30.Pr, 73.43.Lp, 03.65.Vf 



Anyons are two dimensional particles that.unlike bo- 
son or fermions, satisfy exotic statistics [H H, [3, Si- 
These are manifested for abelian anyons by a phase fac- 
tor when they are interchanged or for non-abelian anyons 
by a unitary matrix Q . This complex behavior makes it 
a challenging task to find a representation, mathematical 
or physical, that reproduces these properties. Recently, 
there has been increased interest in these models due 
to their connection to fault-tolerant quantum computa- 
tion 0, [1], and their relation to new states of topolog- 
ically ordered matter Q. They are also of interest in 
the study of multipartite entanglement, due to the so- 
called topological entanglement [l3| required to realize 
anyons on phy sical s y stems. Several proposals have been 
made 11, l3, 13, [l3| for physical systems with anyonic 
behavior. Of particular interest are lattice models, where 
qubits are placed on a two dimensional surface with their 
states representing the vacuum or anyonic populations 
0, [3, [MS [Hi- While for the quantum double models 



we have a simple spin representation of the corresponding 
states [IHIlSiS, 

it has been rather hard to identify 
the states of other models, such as the so-called Ising 
anyonic model. 



Recently, Kitaev [15[ presented a spin lattice Hamilto- 
nian that for different coupling re gim es exhibits abelian 



or non-abelian anyonic behavior [21| , |22| |. The former 



corresponds to the well studied toric code model, while 
the latter corresponds to the Ising anyonic model, whose 
properties have proven difficult to demonstrate. It has 
been shown in the context of the fractional quantum Hall 
effect that states of the Ising non-abelian anyons can be 
built up from states with abelian statistics [23l [2^ . In 
this letter we demonstrate that the topologically entan- 
gled states of a spin lattice capable of supporting the 
anyons of the toric code are sufRcient also to support the 
anyons of the Ising model. This is done by reproducing 
the states of the Ising model through quantum superpo- 
sitions of the toric code states. Our aim is to present the 
mechanism that is responsible for their exotic properties 
without invoking the Hamiltonian description. In this 



scheme the Ising model anyons are located at the end- 
points of strings of operators that belong in the Clifford 
group. Thus, for the first time we give a lattice repre- 
sentation of a non-abelian model that is not a quantum 
double. We achieve this by demonstrating that complex 
topological models can be constructed from simpler mod- 
els, providing a general methodology to perform other 
such mappings. 

To define an anyonic model one needs to give a set of 
possible particle types and their fusion rules. For abelian 
models these rules take the form a x b = c, with only 
one outcome possible for each fusion. For non-abelian 
models, however, the rules take the more general form 
axb = J2c ^ab^' '^ith multiple possible fusion outcomes. 
These rules also determine the quantum dimension, da, of 
each particle, a, with the total quantum dimension of the 
model defined by — X^a^a- Further, one needs the 
R- and _F-matrices that describe the action of a counter 
clockwise exchange of particles and changing the fusion 
order, respectively (see Fig. [Ha,b)). These must satisfy 
the so-called pentagon and hexagon relations which re- 
strict to a finite set of consistent theories [2^ . These con- 
cepts can be summarized using the topological S'-matrix, 
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which can be interpreted as the vacuum-to-vacuum pro- 
cess depicted in Fig. (He). 

Toric code and Ising models:- The toric code model 
consists of four different particle types, the vacuum, 1, 
the anyons e and m and the fermion e. The non-trivial 
fusion rules for these particles are 

exe — mxm — exe ^ 1, exm = e, exe^ni, mxe~e. 

The quantum dimensions are = 1 for all a, hence 
— 4. The R matrices of interest to us are given by 

Rl<^ = {Rl^f = -1, = i?™™ = 1, 
and all of the F-matrices are equal to the identity. 
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FIG. 1: The (a) F matrices, (b) braidings R and (c) S-matrix 
of a general anyonic model. The vertical axis represents time 
running downwards. 



It is possible to represent these anyons on a honey- 
comb lattice with qubits placed on the vertices. Each 
plaquette P is split into two subplaquettes, labeled s to 
the left and p to the right, in order to facilitate the im- 
plementation of the two different species of anyons. We 
identify with the vacuum a state satisfying A.; | ^) = | ^) 
and Bp\^) = | ^) for all s and p, where Ag — erf erf fif erf 
and Bp = aia2<7^ij\ are products of Pauli matrices acting 
on the four qubits of each subplaquette. States having 
any other pattern of eigenvalues can then be identified 
with anyonic populations. For example, an e anyon pair 
is given by the state, | e, e) = erf | ^), corresponding to a 
— 1 eigenvalue of the on the s subplaquettes neighbor- 
ing site i. A string of cr^'s corresponds to two e anyons 
positioned at its endpoints. Similarly the endpoints of 
cr^ and icr^ strings are m and e anyons, respectively. All 
the properties of these anyons can easily be reproduced 
from the representation of their corresponding strings by 
the Pauli operators [2^. 

On the other hand, the Ising model consists of three 
different particle types, the vacuum, 1, the non-abelian 
anyon a and a fermion, tjj. The non-trivial fusion rules 
are 

crXcr = l + '0, T/iX-i/j^l, (T y, ijj = a. 

Again all particles are their own antiparticles. The quan- 
tum dimensions are given by c?i = = 1 and = \f2, 
which also leads to LP' = 4. The R matrices of interest 
to us are 



-1 



This model has the non-trivial F matrix 



t/4 



(1) 



in the basis 1, V'- This shows that a change in the fusion 
order of a particles leads to a superposition of different 
fusion outcomes. 



Superposition Principle:- There is a simple argument 
that shows a relation between the toric code and the Ising 
models. Consider their corresponding S'-matrices, given 
by 
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defined in the basis, 1, e, to, e and 1, cr, ip, respectively. 
Let us consider the equal superpositions of the e and 
TO particle loops of the toric code. Then the following 
relations hold 
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Further, note that = d^ + d^ and that the total quan- 
tum dimensions of the two models are equal [27 1 . 

These observations motivate us to identify a a particle 
loop of the Ising anyon model with these superposed e 
and TO loops, and to identify the fermions of the two mod- 
els. Using this superposition principle we can demon- 
strate fusion, braiding and statistical characteristics of 
the Ising anyon model. In addition, we employ an ancil- 
lary system that provides the correct chirality, and will 
at times require its state to be an entangled state with 
the lattice. 

The state of a superposition of an e string and an m 
string which end in the same plaquettes, may be written 
as 



|o-i,o-2;j; 



1 
71 



(I 61,62) + j |toi,TO2)), 



(3) 



where, for example, we have used | 61,62) to denote the 
state of an e string with one endpoint a plaquette la- 
beled P = 1 and the other in P = 2 (Fig. This can 
be viewed as a new string whose endpoints can reproduce 
the behavior of the a anyons, and so have been labelled 
as such. The relative sign j is a non-local property of the 
string which cannot be determined by local observations 
of the endpoints. It can be changed by braiding oper- 
ations which may act locally around endpoints, and so 
is not topologically protected. The fermions of the toric 
code will reproduce the behaviour of those of the Ising 
model, so we will identify e strings with -0 strings. 
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The movement of a cr must be performed in such a way 
that it does not affect the superposition that encodes the 
non-abehan character of the anyons (2^. This can be 
done by using a qubit ancilla, initially in state | 0)^, and 
the controlled operations 

Cs = + a, + i(l-^)®cr^, 

A = fTr®|0)(0|, + af (4) 

Applying Cs entangles the s plaquette at the endpoint of 
a a string with the ancilla. The operation Di may then 
be applied between the ancilla and the lattice qubit i to 
extend the string one step. To unentangle the ancilla Cg' 
is applied, using the plaquette s' at the new endpoint 
of the (7 string. This method of extending the strings 
is local, allowing the interpretation of their endpoints as 
particles. Also, these operators reproduce the braiding 
statistics of the toric code. 




FIG. 2: The state of a cr string with endpoints in two plaque- 
ttes of the honeycomb lattice can be described by a superpo- 
sition of e and m strings. The relative ± sign is a non-local 
property that cannot be accessed by measurements at either 
endpoint. 

Because we represent particles as endpoints of strings, 
we do not consider any process which cannot be described 
purely in terms of them. So we restrict ourselves to only 
considering the fusion processes in the Ising model that 
can be thought of as two strings fusing to form another. 
This process will be referred to as the fusion of strings. 
We further restrict that the composite string belongs to 
the vacuum sector. As an example, the state of two a 
strings can be written as 

I (0-l,cr2;j)(CT3,CT4;fc)) = 

2 (I ei, 62, 63, 64) -I- j/c I mi, 7712,7713, TO4) + 

k I 61,62,7713,7714) + j 1 7111,7712,63,64)), (5) 

where the relative signs are given by j^k € {— 
The fusion of the strings is achieved by the fusion of the 
particles residing at plaquettes 1 and 3 and of those at 
2 and 4 [2^. The endpoints in the first two terms will 
each behave as the vacuum since each composite object 
is made up of either two 6's or two m's, and similarly the 
second two terms will each give a fermion string. If j ~ k 
the result of the fusion is then 

I (o-l,cr2;j)(o-3,CT4;i)) = 4^(1 ll,3,l2,4> + J I V'l,3V'2,4) ) , 



where we have used | li,3, 12,4) = (161,62,63,64) -I- 
I mi, 7112, 7113, 7ii4))/\/2 to denote the terms that fuse 
to vacuum and \4'i,34'2a) = (| 61, 62, ms, m4) + 
I mi, m2, 63, e4))/\/2 to denote the terms that fuse to a 
fermion string. 

Let us identify | cr, a; ±) with pairs belonging to the 
vacuum and fermion sectors, respectively. The above re- 
sult then reproduces the F matrix ([T]) , and therefore the 
fusion properties, of the Ising anyon model. The fusion 
of a cr string belonging to the vacuum sector with one 
belonging to the fermion sector does not result in a com- 
posite string, and so we need not consider cases where 

By considering the decompositions of the a and ip par- 
ticles in terms of the toric code particles we can show 
that they satisfy the Ising model braiding rules. For ex- 
ample, let us consider the exchange of two ij^'s. Since 
these are identified with the e's of the toric code they 
will have the same fermionic behavior. Also, since the 
braiding of an 6 or an m around an e results in a phase 
factor of —1, so does the braiding of a cr around a ip. 
Let us also consider the braiding of two of the a parti- 
cles, such as those in Eq. ([5]). Braiding the a residing at 
plaquette 1 around that at 3 results in a change of the 
relative sign for both cr strings, and so a change also of 
the relative sign between the vacuum and fermion strings 
in the fusion outcome. From this we infer the R matrices 
(i?J'^)2 = 1 and {R^")^ = -1. These are similar to those 
of the Ising model, except that a complex phase factor is 
missing. This required phase differs for counter clockwise 
and clockwise braidings, e~^'^^'^ for the former and 6*^^/^ 
for the latter. Since R — R'' for the toric code particles, 
the lattice does not distinguish between counter clockwise 
and clockwise evolutions. A framing is therefore pro- 
posed for the cr particles to make this distinction and to 
encode the chirality on an ancillary system. 

We allocate two framings to each cr particle, one to 
the left (l) and one to the right (r). Each of them has an 
ancillary qubit, initially in the zero state, | 0); | 0)^. When 
the particle moves the framings move with it, performing 
the operation 

i?, = l,®|+)(+|+iaf®|-)(-|. (7) 

between their ancillary qubits and the lattice sites, i, to 
the left and right of the particle. This creates super- 
positions of the vacuum and a fermion on the lattice, 
controlled on the ancilla state. When the loops are com- 
plete the framings act trivially on the lattice, but may 
cause a bit flip on the ancilla depending on whether the 
fermion loop acquired a —1 by crossing a a string. Af- 
ter each loop the ancillary qubits are measured and the 
operations e^'"l^a% and e~^'^l^af applied for the results 
I 0); 1 1)^ and I l)j I 0)^ respectively. These assign a phase 
and reset the qubits. The state |0); 1 1)^, for example, 
is assigned e^'"l^ since it is the result of either a counter 
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clockwise loop that encloses no other a particle or a clock- 
wise loop which does enclose a a particle. In the former 
case this phase comes from the fact that the loop causes 
the extended object of the tr particle and framing to un- 
dergo a counterclockwise twist of 27r. This must therefore 
be assigned the phase e"/*, due to a topological spin. In 
the latter case the phase comes from both a clockwise 
braiding and a twist, e*'^/^e~*'^/® = e™!'^ . The consis- 
tency of this framing can be verified in Fig. [3l where a 
complete set of elementary cases have been considered. 

The phase factor required for the i? matrix is that for 
a braiding in which a <j particle performs a loop around 
another particle without twisting. So the twists must 
be removed from the above loops in order to obtain the 
corresponding evolutions. This can be done by following 
all loops with a twist alone in the opposite direction. By 
this two stage process the framing applies the required 
phase of e^*'^/"' for a counter clockwise braiding. When 
the phase is inserted it gives the i? matrix required for 
the consistency of the Ising model. 




(c) 



FIG. 3: (Color online) The four possible loops for a a particle 
that start and finish at the marked points, where the framing 
is depicted. The loops to the top are clockwise and those 
to the bottom are counter clockwise . The loops to the left 
enclose no other a particle and those to the right do. In (a) 
and (d) the left framing crosses a a string once and the right 
framing does not cross or crosses twice, resulting in a bit flip 
on the ancillary qubit for the left framing only. In (b) and (c) 
the situation is reversed. 

We may define the plaquette operator Wp — AgBp on 
the plaquettes of the honeycomb lattice, where s and p 
are the subplaquettes of P. This detects whether an e or 
an m is present on P without distinguishing between the 
two. It can therefore detect the presence of a a particle 
as defined above, without collapsing or otherwise chang- 
ing the superposition. The Wp's defined in this way are 
equivalent to those of Kitaev's honeycomb lattice model 
hj. This means that e and m particles and their super- 
positions correspond to the vortex-like excitations of this 
model, as one would expect for the Ising tr's. Similarly 
the fermions, which are not detected by Wp, correspond 
to non- vortex like excitations. 

Conclusions:- In this letter we demonstrated that the 
superposition of the states of the toric code, together with 



the appropriate framing, can reproduce the fusion, braid- 
ing and statistical properties of the Ising model. This is 
a surprising connection between an abelian and a non- 
abelian anyonic model that reveals the non-local char- 
acter of the latter. It also gives the first lattice repre- 
sentation of a non-abelian model that is not a quantum 
double. It is an exciting possibility to verify if such a 
relation holds between other models, to derive a Hamil- 
tonian that has these states at its low energy spectrum or 
to implement these non-abelian states in the laboratory 

[13, iH. 
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